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We have recorded extreme ultraviolet spectra from to ions using a new flat 

field spectrometer installed at the Shanghai high temperature superconducting electron beam 
ion trap. The spectra were recorded at beam energies ranging between 200 eV and 400 eV 
and showed spectral lines/transition arrays in the 170 - 260 A region. The charge states and 
spectra transitions were identified by comparison with calculations using a detailed relativis¬ 
tic configuration interaction method and collisional-radiative model, both incorporated in 
the Flexible Atomic Code. Atomic structure calculations showed that the dominant emis¬ 
sion arises from 5d —>■ 5p and 5p —>■ 5s transitions. The work also identihed the ground-state 
configuration of as 4/^^5s^ both theoretically and experimentally. 
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I. INTRODUCTION 


There has been a strong interest in tungsten spectroscopy due to the potential use of this 
element as a plasma facing material in the the International Thermonuclear Experimental Reactor 
(ITER) tokamak 1, especially in the divertor region Q]. Spectra diagnostics of the ITER 
divertor region will therefore require a large amount of, as yet, unavailable tungsten atomic data. 
The ITER divertor soft x-ray spectrometer will operate in the 150 - 400 A region 4| where only 
108 lines are known Q]. Of these 102 lines originate from and One line is 

identified as from which is too highly ionized to be observed in the divertor region. The 

plasma temperature in the divertor region of ITER is expected to be in the region of a few to a few 
hundred eV which will result in tungsten in charges states of up to around 28-I-. The remaining 5 
tungsten lines in this wavelength region are reported to be from , but might be mis-identified 
[g]. In contrast to this, the diagnostics of the core plasma region has a stronger support. 


since 


considerably more spectroscopic work, both experimental and theoretical, has been reported for 
the important charge states, (see the review by Kramida ^). 

In previous papers we have investigated visible Ml transitions in Q], 0, W26+ Q, 

W27+ jq| an extension to the silver iso-electronic sequence [3,Q) and W28+ Q . In the 

present work we report on a study of tungsten ions in the wavelength region of interest to ITER 
divertor diagnostics, i.e. 200 - 400 A. Using the Shanghai high temperature super conducting Elec¬ 
tron Beam Ion Trap (SH-HtscEBIT), which was designed for low electron beam energy operation, 
we observed spectra from through to The spectra were recorded using a recently 

developed high resolution flat field spectrometer [l^ . From calculations of the relevant atomic 
structure and also collisional radiative model, using the FAC code [l^, we have identified lines and 
spectral features as originating from hd —> 5p and 5p ^ 5s transitions. We will also discuss the 
interesting case of in more detail, due to the controversy over the classification of the ground 

state which has existed since the calculations by Curtis and Ellis in 1980 Q- 


II. EXPERIMENTAL METHOD 


. This 

EBIT is capable of operating in the range of electron beam energies between 30 and 4000 eV. The 
magnetic field, which is created by liquid nitrogen temperature superconducting coils, compresses 
the beam radius to 60 /j,m. The background vacuum pressure in the trap center is estimated 


The experiment was carried out using the SH-HtscEBIT, which is described elsewhere 
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TABLE 1. The calibration lines used in this work, from the background oxygen(0 IV 
nitrogen(N IV 25|, N V [^) lines. 


231, O V [2J) and 


Spectra OV OV NV NIV OIV 

Transition iPi - ^Sq ^Di.s - ^ ^ 

Wavelength(A) 172.169 192.799 209.303 225.210 238.570 


to be lower than torr, which makes it possible to produce tungsten ions mainly through 

electron collisional ionization with negligible influence from charge exchange. The spectra were 
recorded by utilizing a high resolution grazing-incidence flat field spectrometer, which covers the 
range of 10 to 500 A and reaches a resolving power of above 800 151. In order to eliminate light 


from the hot cathode, a 4500 A thick aluminum foil was used and mounted on the window of the 
mini ultra high vacuum gate valve between the SH-HtscEBIT and the spectrometer. Due to the 
aluminum L absorption edge, only wavelengths longer than 171 A can pass through this foil. For 
the present experiment, a Shimazdu varied-line-spacing (VLS) grating (1200 1/mm, part number: 
001-0659) 19 -211 was used and an Andor CCD camera (model number: DO936N-00W-#BN) was 
placed at different positions to record different wavelength regions. Tungsten ions were obtained by 
injecting W(CO)6, a volatile compound with a high vapour pressure at room temperature 2^. The 
experiments were done using electron beam energies ranging from 200 to 400 eV, in steps of about 
20 eV, and a beam current of 8.1 mA. The spectrometer was calibrated by several background 
oxygen and nitrogen lines (see table H]) . 


III. DESCRIPTION OF THE CALCULATION 


In order to support the identification of the recorded spectra and to predict the charge state 
distributions we performed calculations using the Flexible Atomic Code, FAC vl.1.1. [l^. This is 


an integrated software package 


27h 29I| producing both structure and scattering data, including e.g. 


energy levels, radiative transition rates, collisional excitation and ionization by electron impact. 
FAC also includes a collisional radiative(CR) model to produce synthetic spectra for plasmas under 
different physical conditions. The atomic structure calculation in FAC is based on a relativistic 
configuration interaction (RCI) model using independent particle basis wavefunctions, derived from 
a local central potential. The orbitals are optimized in a self-consistent-field (SCF) iterative pro¬ 
cedure. Relativistic effects are included through the Dirac-Coulomb Hamiltonian being used for 
the optimization. Corrections from Breit interacting and higher order QED effects, e.g.vacuum 
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TABLE II. Configurations included in the CR-model for to (see text), together with the resulting 

number of excited levels. We also give the predicted ground state of each charge state in LS coupling and 

jj-coupling. 

ions: + 

W33+ 

w33 + 

W34+ 

W33+ 

# levels: 2538 

1769 

1196 

3896 

3655 

4f^‘'5s^Sp 

4/1*5s3 

4/l*5s 

4/34 

4/33 

4f^^5s5p^ 

4/3*5s5p 

4/3^5p 

4/3358 

4/3358 

4/l'‘5sSp5d 

4/l*5p3 

4/1^5d 

4/335P 

4/335P 

4/1^5s^5p^ 

4/'-*5s5d 

4/335s2 

4/335d 

4/335d 


4/^^5p5(i 

4/335s5p 

4/335s3 

4 / 335*3 

4/135sSp3 

4f^^5d^ 

4/335p2 

4/335s5p 

4/33585p 

4f^^5s5p^5d 

4/335s25p 

4/335s5d 

4/335p2 

4/335p3 

4/135a25p3 

4/335s25^ 

4/325s25p 

4/335s5d 

4/^^ 5s5«i 

4f^^5B^5p‘^Sd 

4/335s5p2 

4 / 335 ^ 25 ^ 

4/335s35p 

4f^°5s^5p 


4/335s5ci2 

4/335s5p3 

4/335s35d 

4f^°5s^5d 


4/^^5s5p5<i 

4/335825p2 

4/335s5p3 



4/'-35s25p2 


4/305s35p3 



4/^^5s^5p5d 





4/335825^2 





ground 

4/335835p3 

4/3*5s 3 3gg 

4jrl3i-„2 2 ^ 

4/ 5s F 7/2 

4/33683 3He 

• /IIc„2 4 t 

4/ os I 15/2 

state: 

[(4/^^) 7/2 ^ 2 ) 0^772 

[4/3^5s3]o 

[(4/^^)7/2 5s^]7/2 

[[(4/®;2)o(4/®/2)6l65*^]6 

[[{4/5/2)o(4/7/2)i5/2]i5/2 5s^]i5/2 


polarization and self-energy, are included in a final Cl-calculation. This approach yields en ergy 


levels, radiative and autoionization rates, and collision strength from the final wavefunctions 


M- 


The CR-model has been shown to be a reliable tool for the analysis of optically thin plasmas 


Q- 


such as those in an EBIT [SOU. CR-model can provide information on level populations and spec¬ 
tral line intensities and therefore be used for the generation of synthetic spectra. In this model, 
we have considered electron-impact excitation and de-excitation, together with radiative transi¬ 
tion, but neglected other effects, e.g. charge exchange and radiative recombination, since they are 
less important in dilute EBIT plasma. For a given excited level, with the normalization condition 
= 1, the population was obtained by solving quasi-stationary-state rate equations = 0, 

where 


j>i j<i j>i j<i j>i j<i 

(1) 

In this equation and Nj are the radiative transition rate from level j to i, the 

collisional excitation rate from j to i, the collisional de-excitation rate from j to i and the population 
of level j, respectively. 

Generally the accuracy of our calculations depends on the size and completeness of the atomic 
model used, which is determined by the choice of included configurations. Given the complexity 
of tungsten ions in the range of charge state of interest here we were forced to consider a limited 
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TABLE III. The ionization potential of to 31|. 


charge state 







ionization energy (eV) 

208.9 ± 1.2 

231.6 ±1.2 

258.2 ± 1.2 

290.7 ± 1.2 

325.3 ± 1.5 

361.9 ±1.5 


number of levels in our CR-model. One important consideration was to include all important 


metastable levels since they have 


ong lifetimes, which in turn also gives high population, even 


14i |. These properties of the metastable levels can lead to two 


relative to the ground state 
step excitation and ionisation 1^. In to we will deal with configurations with open 

4/-subshells, which are known to give rise to several metastable levels due to the high angular 
momenta involved. To give models of reasonable complexities, we restricted our model to include 
configurations created by single and double excitation from 4f, 5s and 5p subshells, as shown in 
table HU where we also give the total number of states included for each ion. The synthetic spectra 
were create using CR-models for beam energies of 230, 270, 310, 350 and 390 eV, respectively, 
for to at an electron density of 5 x 10^® cm“^, which is a typical electron density for 

the SH-HtscEBIT. The calculated lines are convolved with a Gaussian line shape to match the 
spectrometer resolution of 850. 


IV. RESULTS AND DISCUSSIONS 


The measurements reported in this work scanned the wavelength region of 100 - 400 A, but all 
identified lines fall in the 160 - 270 A region. The spectral development as a function of electron 
beam energy, and therefore the charge state of the tungsten-ions, is shown in figure [TJ 

Based on the ionization potentials given in table HU we assign the features labelled 1, 2, 3 and 
4, which appear at electron beam energies between 200 and 230 eV, to the spectrum of In 

the same way, we conclude that the features labelled 5, 6 and 7 belong to lines, the transition 

array 8 is from , 9 is from and 10 is from The identification of the arrays 

are further supported by the fact that they appear at the same beam energy as the visible lines 
in the same ion, as identified in 8|. In figure [2] we compare these experimental spectra with our 
synthetic ones from the CR-model spectra at some selected electron energies, where a single charge 
state is dominating for each case. We also list our computed ground states of to which 

agrees wth earlier results by Kramida and Shirai 311. 


The synthetic spectra give us important guidance when identifying the experimental results, but 
there are two clear and systematic deviations. First, there is a shift towards longer wavelengths 
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FIG. 1. Spectra taken at SH-HtscEBIT at 11 electron beam energies from 200 to 400 eV, with a beam 
current of 8.1 mA, while injecting W(CO)6. Based on the ionization potentials given in table Hill the lines 
labeled 1, 2, 3 and 4 belong to lines, 5, 6 and 7 belong to lines. Transition array 8 is from 

Wi3+, 9 is from and 10 is from Wi^+. 


in the calculations, and second the transition arrays are too wide in the synthetic spectra. It 
may be possible to correct for these differences by including more correlation, such as core-core 
(CC) correlation 111], but since we are aiming for CR-modelling with a large number of atomic 
properties, we have to limit the size of our atomic structure model. It is clear, though, that the 
features of the synthetic spectra can be used to identify the transition arrays appearing for each of 
the tungsten charge states. 

The lines are for all ions concentrated in the two wavelength regions (see figure[2|) 160 - 190 A and 
210 - 260 A. The former corresponds to 5d — 5p transitions, while the latter arises from 5p —)• 5s 
arrays. The results are in agreement with Suzuki’s work 3^, in which their calculation showed 


5p —)• 5s transitions dominating the spectrum around 250 A for the tungsten charge states of 
Our detailed identihcations are listed in table nYi and we will discuss them for each 
ion in detail in the following paragraphs. 
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160 170 180 190 200 210 220 230 240 250 260 270 

Wavelength (A) 


FIG. 2. Experimental (exp) and synthetic (cal) EBIT-spectra of to at electron energies of 230, 

270, 310, 350 and 390 eV, respectively. The calculated lines were convoluted with a line width to resemble 
an instrument resolution of 850. The numbers refer to the same transitions as in figure [1] and table IIV| Lines 
labeled “a” and “b” in the synthetic spectra are not possible to observe in the experiment due to the A1 foil 
transmission properties (giving rise to the cut-off around 170 A). 
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TABLE IV. Identified lines and arrays in The labels in the second column refer to figure [H 

The wavelengths are in unit of A. The uncertainties of the experimental wavelengths in transition arrays 


are around 0.15 A and mainly due to the statistical uncertainty of the wavelength calibration and fitting. 


ion label experimental 

theoretical 

transitions 

type 

wavelength 

wavelength 




W“+ 










1 175.6 - 180.4 183.3 - 189.0 4/i®5s^5p5d ->■ 4:f^5s^5p'^ array 

2 246.23 ± 0.12 250.84 [[(4/i®)s/25s]3(5p?/25p3/2)3/2]3/2 ^ [[(4/^")5/25s=]5/2(5p?/2)o]5/2 line 

3 252.25 ± 0.18 257.46 [[(4/i®)s/25s]2(5p?/25p3/2)3/2]7/2 ^ [[(4/^")5/25s2]5/2(5p?/2)o]5/2 line 

4 254.09 ± 0.17 261.98 [[(4/^")5/25s]3(5p;/25p3/2)3/2]5/2 ^ [[(4/”)5/25s^]5/2(5p;/2)o]5/2 line 

5 174.8 - 176.8 173.8 - 178.0 4f^5s^5p5d Af^5s'^5p'^ array 

6 243.27 ± 0.13 243.29 mf^)7/25s]45pi/2]9/25p3/2]4 ^ [(4/i®)7/25s=5pi/2]4 line 

7 249.19 ± 0.13 257.43 _ [[(4/^")7/25s]35pi/2]7/25p3/2]5 ^ [(4/^^)7/25s"5pi/2]4 line 

a 173.3 - 178.08 4f^5s^5d 4f^5s^5p array 

f 4/13 5s5p4/1® 5s^ 

8 232.7 - 240.9 225.8 - 245.9 < array 

[ 4/12 5s5p^ 4/i2 5s2 5p 

b 166.6 - 171.3 4/ii5s25d 4/ii5s25p array 

f 4/12 5s5p4/12 5^2 

9 225.4 - 233.7 221.5 - 242.7 < array 

[ 4/ii5s5p2 4/ii5s25p 

10 214.8 - 226.9 218.4 - 234.5 4/ii5s5p ^ 4/ii5s2 array 


a. Our predicted ground state of is 4 /^ 25525^2 ([[( 4 /|/ 2 )o( 4 /J/ 2 ) 7 / 2 ] 7 / 25 s^( 5 p^/ 2 )o] 7/2 
in jj-coupling). The experimental and synthetic spectra are shown in figure[2]at E = 230eV. Accord¬ 
ing to the spectra, four transition arrays or lines are determined to be from this ion. 4 /^ 25525 ^ 5 ^^ 

^ 4/^25525^2 transitions contribute to the 175.6 - 180.4 A array(labeled as “1” in figure [2]) and 
the three longer-wavelength lines we identify as 4/^2555^3 ^^-135525^2 transitions. 

b. In the present work the ground state of is [4/^‘^5s^]o. We also compared our 


computational results with available data in 
various results from 


33|. The differences between our RCI results and the 


33|, i.e. RMBPTl(First-order relativistic many-body perturbation theory), 


RMBPT2(Second-order relativistic many-body perturbation theory) and COWAN(Hartree-Fock 
relativistic method) values of excitation energies are 1-5%, 0-3% and 1-4%, respectively. Moreover, 
our results are in better agreement with the RMBPT2 results and most of the excitation energies 
agree within 2%. We would like to point out that in spite of the fact that the ground state s ugg ested 


from our calculations is the same as the ones given by Kramida and Safronova’s works [31|, l33l | . 
this is only a tentative identification. The level splitting of the 4/^^5s^5p^ configuration is very 
sensitive to the computational model, and the inclusion of more correlation migth push the lowest 
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FIG. 3. Energy level structure of “a” refers to energy levels of the and 4/^^5s5p^ 

configurations, while “b” gives energy levels of the Af^^bsbd and Lines of observed in this 

work are mainly from Af^'^bs^bd —>■ 4/^^5s^5p, 4/^^5s5p ^ 4/^^5s^ and 4/^^5s5p^ —>■ 4/^^5s^5p transitions. 
The “Curtis and Ellis” lines (see text) are marked in the right portion of the figure. 


level in this configuration below the 4/^^5s^ ^Sq level. A larger calculations will be needed to 
resolve this issue. The experimental and synthetic spectra are shown in figure [2] for an energy of 
E = 270eV. The transition array “5” in the wavelength region of 174.8 - 176.8 A in the experimental 
spectra(see figure [2|) is identihed to belong to the Af^‘^5s‘^5p5d 4/^^5s^5p^ array, while the other 

two lines at 243.3 A and 249.2 A are bp —>■ 5s transitions, outside a 4/y^2"Subshell. 

c. The experimental and synthetic spectra of this Promethium-like ion are shown in 

figure [2] at E = 310eV. There is still some remaining controversy concerning the ground-state of 
because of the effect of the “collapse” of the 4/-orbital, which will lead to a change in 
ground-state configuration of promethium-like ions from 4/^^5s^ to 4/^^5s, leading to a simpler 
semi-one-electron spectrum, with two strong resonance lines between the 5s and bp levels (see 
figure [3]). In 1980, Curtis and Ellis predicted that this change would happen for Z < 74, by using 


the Hartree-Fork method 


Q- 


Later, in 1983, Theodosiou indicated that the change occurred first 


at Z = 78 by using the Dirac-Fork method 


n 


34l | . Recent calculations |3,5l437l| support the prediction 
of Theodosiou. In 2003, Hutton et al. reported a detection of the 5p ^ 5s resonance lines in 
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TABLE V. Wavelength for the ^Pi/ 2 , 3/2 4/^^5s ^Si /2 transitions in The values are in unit 


transitions 

RCI 

RMBPT^ 

RMBETB*" 

COWAN'= 

MR-MP'i 

exp® 

^P 1/2 ^Si /2 

367.52 

370.92 

371.04 

367.95 

369.14 

365.3 

^P3/2 ^Si/2 

257.83 

260.94 

261.12 

260.83 

259.62 

258.2 


Second-order approximation with a Dirac-Fock (DF) potentials 
Second-order approximation with a Breit-Dirac-Fock (BDF) potentia 
Hartree-Fock relativistic method implemented in tlp-COWAN code | 
Multireference M(^ller-Plesset (MR-MP^ approach [3g. 

® Experimental wavelengths from Ref. |38l] . 


3 - 


a spectrum from the Berlin-EBIT 


38l |. which would support the model by Curtis and Ellis, and 


contradicting later theoretical predictions. However, the performance of this high-energy EBIT, 
working at only a few hundreds of eV energy, is quite debatable and leads to an uncertain and 
tentative identification. 


The wavelengths of the "^Pi /2 —t 4/^^5s "^Si /2 and 4/^^5p "^P 3/2 —t 4/^^5s "^Si /2 

transitions for are listed in table In the present work our calculations predict the 

4/^^5p ^Pi/ 2 , 3/2 4/^^5s ^Si /2 transition lines at 367.52 A and 257.83 A. The biggest difference 

between our RCI and other theoretical results for the ^Pl/2 ^Si/2 and ^P 3/2 —1- ^Si /2 transi¬ 

tions are 0.96% and 1.27%, respectively. The difference between our RCI and the experimental 
results for the ^Pi /2 —>• ^Si /2 and ^P 3/2 —)• ^Si /2 transitions are 0.14% and 0.60%, respectively. 
The ground state configuration is predicted to be 4/^^5s^ and our CR-model do not indicate any 
strong Af^^bp —)■ 4/^^5s lines. This is in agreement with our experiment where we did not observe 
any strong resonance line around the predicted position (see table ED with a simple structure. 
Figure [3] shows the level structure of Clearly there are numerous levels lower in energy than 

the 4/^^5s and 4/^^5p. This leads to a low population of t 
as opposed to the strong lines predicted by Curtis and Ellis 
and where the simplest configurations, with closed 4/^^-subshell configurations, are not the 

ground state. The occurrence of simple spectra, with few resonance lines, will therefore occur for 
higher Z in all the corresponding iso-electronic sequences. To further support the identification 
of 4/^^5s^ as the ground-state configuration of we recently observed visible MI transitions 

within the 4/^^5s^ and 4/^^5s^5p configurations using the SH-HtscEBIT The transition array 
“8” in figure [ 2 ] observed at around 232.7 - 240.9 A was identified as bp —)■ 5s arrays outside a 4/^^ 
or 4/^^ subshell (see table lYD, which coincide with the results of Safronova [^. A second strong 


le 4f^^5» levels and hence weak lines 

□ 

I7l |. This situation is similar to 












11 


transition array “a” in figure [H positioned at 173.3 - 178.08 A arises from Sd —>• 5p transitions, 
outside a 4/^^ subshell, but is unobservable in the experimental spectrum, due to the wavelength 
cut-off of the Aluminium foil. 

d. We predict the ground state of to be ([[(4/|/2)o(4/7/2)6]65s^]6 in jj- 

coupling). The experimental and synthetic spectra are shown in figure [2] at E = 350eV. The 
transition array “9” in figure [2] covering the wavelength region of 225.4 - 233.7 A is suggested to be 
from 5p 5s transitions outside a 4/^^ or 4/^^ subshell (see table lYl). The shorter wavelengths 
“b” in hgure[2]at around the 166.6 - 171.3 A wavelength region arises from 4/^^5s^5(i 4f^^5s^5p 

transitions. Again these lines did not appear in the experimental spectrum due to the A1 foil 
transmission properties. 

e. Our predicted ground state of is ([[(4/|/2)o(4/7/2)i5/2]i5/25s^]i5/2 in 

jj-coupling). The experimental and synthetic spectra are shown in figure [2] at E = 390eV. The 
transition array “10” at 214.8 - 226.9 A in figure [2] are identified as from 4/^^5s5p —>■ 4/^^5s^. 


V. CONCLUSION 


To conclude, we have studied EUV spectra from to ions experimentally and the¬ 

oretically. The experiments were carried out using a new flat held spectrometer installed at the 
SH-HtscEBIT. The spectra were recorded at electron beam energies ranging between 200 and 400 
eV and spectral features were recorded in the 170 - 260 A region. In order to identify the charge 
states and spectral features we performed detailed RCI calculations of the relevant atomic struc¬ 
ture. The data was then used in a CR-model of the EBIT. Both calculations are a part of the 
EAC code. Although there were systematic shifts in the calculated wavelengths and the widths of 
the spectral features were a systematically wider than the corresponding experimental arrays, the 
synthetic spectra are in good enough overall agreement with experiment, to facilitate a number 
of identifications of arrays and lines. The spectral features in the 160 - 190 A region are from 
5d —>■ 5p transitions whereas 5p —>■ 5s falls in the 210 - 260 A region, for all the ions considered. 
This work also identihed the ground-state conhguration of to be 4/^^5s^, supported by both 

experimental spectra and theoretical calculations. 
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